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ABSTRACT 

Far-ultraviolet (FUV) radiation plays an important role in determining chemical abundances in pro- 
toplanetary disks. H I Lyman a is suspected to be the dominant component of the FUV emission from 
Classical T Tauri Stars (CTTSs), but is difficult to measure directly due to circumstellar and interstel- 
lar H I absorption. To better characterize the intrinsic Lyman a radiation, we present FUV spectra 
of 14 CTTSs taken with the Hubble Space Telescope COS and STIS instruments. H2 fluorescence, 
commonly seen in the spectra of CTTSs, is excited by Lyman a photons, providing an indirect mea- 
sure of the Lyman a flux incident upon the warm disk surface. We use observed H2 progression fluxes 
to reconstruct the CTTS Lyman a profiles. The Lyman a flux correlates with total measured FUV 
flux, in agreement with an accretion-related source of FUV emission. With a geometry-independent 
analysis, we confirm that in accreting T Tauri systems Lyman a radiation dominates the FUV flux 
(~1150 A - 1700 A). In the systems surveyed this one line comprises 70 - 90 % of the total FUV flux. 
Subject headings: stars: pre-main sequence 



1. INTRODUCTION 

One of the most influential factors in determining the 
local composition and physical state of the gas in pro- 
toplanetary disks is the ultraviolet (UV) radiation field. 
Strong UV radiation is produced from hot gas in the 
magnetically active atmospheres of the central star, and 
at the magnetospheric shock where disk material collides 



with the stellar atmosphere (Giinther & Schmitt 



Sacco et al. 2008 1 . UV photo ns have a profound inf luence 



2008 



on the gas heating (jJonkheid et al.| |2004 Nomura et al 



20071 i Woitke et al. 2009) an d disk gas chemistry (|Aikawa 
k Herbst|l999||Betliell fe Bergin|2009[ |Fogel et al.|2011[ ). 

Surveys of the U V radiation field have been able to char- 
acterize the general spectral structure of the radiation 



(Valenti et al. 2000 Yang et al. 2012 1. H I Lyman a 



(LyaJ is a significant contributor to the UV radiation, 
comprising as much as 80% of the far-UV (FUV) emis- 
sion p roduced in the stellar atmosphere and a ccretion 
shock ( |Bergin et aL]|2003[ |Herczeg et ~aT1|2004| ). While 
the Lya radiation is intense close to the star, atomic hy- 
drogen in the disk atmosphere at larger radii isotropically 
scatters Lya photons which provides greater penetration 
into the molecular disk. This is quite different from UV 
photons at other wavelengths that are scattered solely 
by dust grains, leading to a Lya dominated radiation 
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field - if L ya ra diation is present ( Herczeg||2005 Bethell 
& Bergin 20111. This is important because Lya disso 



ciates many molecules such as H2O and CH4, and can 
lead to an overabundance of certain species, especially 



CN flBergin et al I |2003t |van Zadelhoff et al.|[2003 



Dishoeck et al.||200fc>p 



van 



While many strong U V transitions such as C IV and 
He II have been studied ( Ingleby et al.||20lTj Yang et ah] 



2012|, an intrinsic Lya profile has only been well char 
acterized in a single Classica l T Tauri Star (CTTS) - 
TW Hya ( |Herczeg et al.|2004j ). Due to circumstellar and 
interstellar H I absorption, it is impossible to directly 
measure the intrinsic Lya emission of most CTTSs. On 
the other hand, H2 emission resulting from Lya photoex- 
citation is prevalent throughout the FUV bandpass, pro- 
viding an indirect probe of the Lya. Several Lyman and 
Werner band transitions reside at wavelengths coincident 
with the Lya profile, and have a relatively simple radia- 
tive transfer. These line s have been identified in previous 
observations of CTTSs (Brown ct al. 198l] |Valenti et aL 
[20001 JArdila et al.|2002[ |Herczeg et al |2^02PTO1 prance 
|et al'poil al). H2 emission therefore provides an oppor- 
tunity to study the circumstellar region where the Lya 
radiation interacts with the disk molecular g well 
as characterize the strength of this key component of the 
FUV radiation field in T Tauri systems. 

The impact of strong Lya radiation on disks has only 
recently be en explored in detailed UV radiation trans- 
fer models (IB ethell & B ergin|2011 ) and chemical models 
( Fogel et al.||2UlTj T Given the influence of these pho- 
tons on gas throughout the disk, quantifying the Lya 
emission in CTTSs is crucial to produce accurate disk 
thermal/chemistry models. In this letter, we reconstruct 
the Lya profile from observed H2 emission lines in the 
FUV spectra of CTTSs. 

2. OBSERVATIONS 

Our sample includes 14 CTTSs, chosen for their strong 
H2 and FUV continuum emission. The targets, as well as 
their distances and extinctions, are listed in Table [j] We 
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adopt targe t and line-of-sight parameters from France | 
et al.| ( |2012| ). Most observations utilized the G130M and 
G160 M medium resolutio n (Av «17 km s _1 ) modes of 



! ((Green et al. 2012|). We also use archival STIS 
3M^ 



COS 

E140M obser vations of TW Hya from the StarCAT spec- 
tral catalog (Ayres 2010[) to compare our results with 
previous analysis (Hcrczcg et al.|[2004 ). 

This sample includes targets from the Taurus-Auriga, 
rj Chamaeleontis, and TW Hya star-forming regions, in 
addition to isolated systems. They are a subset of a full 
CTTS H2 survey presented in more detail in France et al. 
(2012). All the CTTSs in our sample display Lya-excited 
H2 emission lines throughout the COS bandpass with an 
average line width of 46 ± 14 km s _1 . Circumstellar 
and interstellar H I absorbs the Lya profiles to vary- 
ing degrees in each target. Figure [T] shows the observed 
H2 emission and Lya profiles of several targets. An H I 
component absorbs the Lya line core in V4046 Sgr, while 
an outflow absorber attenuates most of the short wave- 
length side of the Lya profile of RU Lupi. Very little 
Lya emission is observed in the DE Tau spectra. How- 
ever, each target shows abund ant Lya excited H2 emis- 
sion lines ( |France et al. 2012), demonstrating that the 
observed Lya profile can vastly underrepresent the local 
Lya radiation field at the disk surface. 

3. ANALYSIS 



Wood et al | ( |2002[ ) and |Herczeg et al.| ( |2004[ ) developed 
a method of Lya reconstruction using photo-excited H2 
emission lines in the UV spectra of Mira B and TW Hya. 
The line fluxes were used in fluorescence models to de- 
termine the intrinsic Lya profiles and H2 properties of 
each object. To make an analysis of a larger sample 
tractable, we employ a similar technique with only the 
brightest H2 progressions. A progression refers to the 
cascade of electronic transitions originating from an in- 
dividual rotational-vibrational level of the excited elec- 
tronic state. Similar to the previous studies, we first 
measure the H 2 lines, and then model the Lya profiles. 

3.1. H2 Emission Spectra and the Absorbed Lya Flux 

We employed a multi-Gaussian IDL line-fitting code, 
optimized for COS emission line spectra, to measure the 
total flux from Lya-pumped H2. This code assumes a 
Gaussian line-shape convolved with a wavelength depen- 
dent line spread function (LSF), then uses the MPFIT 



routine to m inimize \ between the fit and data (Mark 



wardt 20091. An unconvolved Gaussian was used for 
TW Hya observed with STIS. To reconstruct the local 
Lya profile incident upon the molecular disk surface, 
we measured the total flux from 12 fluorescent progres- 
sions excited by Lya. We chose to restrict the emis- 
sion line analysis to the 1395 - 1640 A range to mitigate 
the effects of H2 self-absorption, w hich are strongest at 
A < 1400 A ( jHerczeg et al1|2004| ). The brightest, un- 
blended lines from 12 progressions in the 1395 - 1640 A 
bandpass (up to 38 lines) were fitted. Wc refer the reader 
to France et al. (2012) for additional information on the 
H2 emission characteristics of our sample, including pro- 
gression IDs for all lines measured. 
The total flux from progression m is given by 



F, 



N ^ 



B rl 



where F mn is the reddening corrected, integrated H2 
emission line flux (in units of ergs cm~ 2 s _1 ) from rovi- 

brational state m (= [v , J ]) m the B 1 ^ electronic state 
to n (= [v ,J }) in the ground electronic state, A 1 E+. 
B m n is the branching ratio for each transition in a pro- 
gression, defined as the ratio of the Einstein A-valne 
for a given transition m — > n (A v >j>^ v "j>i) to the to- 
tal transition rate out of state m, including tra nsitions 



to bound states and the vibrational continuum ( Stecher 
I fc Williams|[T967| |Wood fc Karovska|[2004l |France et al 



|2011b I . N is the number of emission lines measured from 
a given progression. 

We take the error to be the standard deviation of the 
individual measurements of i 7 ' m (H2), as this is the rela- 
tive uncertainty in H 2 fluxes in most cases. These stan- 
dard deviations are typically less than 10% for bright 
progressions, and as high as 20 - 30% in weaker progres- 
sions. The primary uncertainty in H2 luminosities relates 
to the extinction correction (see Ay in Table 1), although 
to zeroeth order this would affect all progressions simi- 
larly as they have been chosen from a narrow range of 
wavelengths. Upper limits on the H 2 emission line fluxes 
of undetectable progressions were determined from the 
standard deviation in a ± 50 km s" 1 region surrounding 
the laboratory wavelength of the transition. 

The flux incident on the H2 (/mc(A m ), in units of ergs 
cm~ 2 s _1 A -1 ) at each H 2 absorbing transition wave- 
length A m was then assumed to be the total progression 
flux divided by the equivalent width (W\) of the absorb- 
ing H2 transition: 



Iinc{Kn) — F m /W\ 



(2) 



(3) 



where each [v,J] state of the ground electronic level is 
populated by a single rotational temperature (Tjj 2 ) and 
column density (iV#,). Non-thermal populations can re- 
sult from dust heating and pho to excitation by the strong 
X-ray and UV radia tion fields (N omura et al.|2007 Gorti 
& Hollenbach 2009). However, a degeneracy exists be- 
tween and 1\h-, in our model, discussed further in 



Section |3.2| As a first order approximation to simplify 
this degeneracy we assume a thermal population. To de- 
termine W\, we used a grid of T# 2 in 100 K increments 
from 1000 K to 5000 K, and iV# 2 in 0.1 dex increments 

Models with T Ha < 1000 K 



from 10 16 to 10 



(1) 



produce a negligible population in the v=2 level of the 
ground electronic state and are therefore not considered. 
For each of the 12 transitions, we calculated a set of Ii nc 
for each pair of T# 2 and Nh 2 ■ Each pair of T# 2 and Nh 2 
therefore creates 12 unique values of Ii nc . 

3.2. Lya Profile Reconstruction 

Our reconstruction technique finds the set of outflow- 
absorbed Lya profiles (L a b, in units of ergs cm -2 s _1 
A -1 ) that best fit the variety of Ii nc values. The model 
assumes Lya emission is created near the stellar surface. 
Some amount of H I lies in between the Lya source and 
the H2, absorbing some of the Lya photons before they 
reach the H 2 . For each target, we create a grid of I a b 
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Table 1 

Target Properties and Fluxes 



Target A v A < F Lya > a < F ab > b FWHM N out v out F FUV 
(pc) (kms- 1 )' (10 19 cm- 2 ) c (km s~ 1 ) c 



AA Tau 


0.5 


140 


35.2 


± 


6.3 


16.5 


± 


2.4 


642 


± 


33 


1.11 


± 


0.31 


-143 


± 


16 


3.5 


BP Tau 


0.5 


1 10 


31.5 


± 


5.5 


20.7 


± 


3.0 


613 


± 


54 


0.44 


± 


0.30 


-131 


± 


29 


14.1 


DE Tau 


0.6 


110 


15.4 


± 


4.5 


7.91 


± 


2.8 


623 


± 


75 


0.92 


± 


0.53 


-156 


± 


28 


3.6 


DM Tau 


0.0 


110 


4.6 


± 


0.5 


3.19 


± 


0.4 


912 


± 


35 


0.39 


± 


0.03 


-89 


± 


•1 


1.1 


GM Aur 


0.1 


110 


12.2 


± 


3.0 


6.19 


± 


1.5 


815 


± 


31 


1.16 


± 


0.32 


-129 


± 


17 


2.9 


HN Tau 


0.5 


110 


12.9 


± 


2.6 


7.08 


± 


1.4 


776 


± 


10 


1.15 


± 


0.49 


-153 


± 


26 


2.8 


LkCa 15 


0.6 


110 


17.2 


± 


2.8 


9.54 


± 


2.0 


665 


± 


39 


0.80 


± 


0.40 


-139 


± 


25 


3.0 


RECX 11 


0.0 


97 


5.2 


± 


0.4 


3.07 


± 


0.6 


573 


± 


63 


0.78 


± 


0.48 


-126 


± 


35 


0.8 


RECX 15 


0.0 


97 


11.6 


± 


1.6 


5.73 


± 


1.0 


848 


± 


54 


0.96 


± 


0.38 


-125 


± 


27 


1.0 


RU Lupi 


0.1 


150 


23.4 


± 


5.5 


10.3 


± 


1.7 


707 


± 


59 


3.14 


± 


1.56 


-219 


± 


52 


7.5 


SU Aur 


0.9 


140 


35.4 


± 


8.2 


14.8 


± 


1.9 


644 


± 


59 


2.34 


± 


1.37 


-177 


± 


47 


6.9 


TW Hya 


0.0 


56 


52.5 


± 


9.1 


38.3 


± 


6.6 


865 


± 


35 


0.29 


± 


0.12 


-77 


± 


16 


41.4 


UX Tau 


0.2 


110 


6.2 


± 


0.5 


4.32 


± 


0.4 


606 


± 


12 


0.18 


± 


0.04 


-64 


± 


7 


1.1 


V4046 Sgr 


0.0 


83 


46.0 


± 


5.1 


24.9 


± 


2.2 


627 


± 


38 


0.65 


± 


0.26 


-108 


± 


17 


7.0 



Fluxes in 10 12 ergs cm 2 s 1 

a Intrinsic, unabsorbed model flux 

b H2-incident, circumstellar H I-absorbed model flux 

c Average parameters of intrinsic model Lya profiles 

d Integrated 1160 - 1695 A flux, excluding hya (1207 - 1222 A) and O I (1300 - 1310 . 



profiles, each with its own single Gaussian emission com- 
ponent and single outflow H I absorber: 



lab — {Chya + Ihya)G Tout 



(4) 



Chya is the average continuum flux per Angstrom near 
Lya. Ihya is the accretion/magnetospheric-generated 
Gaussian emission profile, centered at the stellar radial 
velocity and parameterized by amplitude Ig and FWHM 
(in units of ergs cm -2 s _1 A -1 and km s -1 , respectively). 
The optical depth r out is determined with a Voigt profile 
for the H I Lya transition, characterized by an outflow 
velocity v out and column density N out (in units of km s 
and cm -2 , respectively). This outflow absorption com- 
ponent is essential to adequately fitting the I inc values. 

For each I a b model profile, a x 2 is calculated with each 
set of [Tf/^iVfl-J-dependent I inc . We then exclude any 
sets of Tjj 2 and Ng 2 which yield a x 2 greater than 95.4% 
(2a) from the peak of the x 2 probability distribution. 
The H2 progressions for which only upper limits are cal- 
culated provide additional constraints on the Lya profile; 
we reject any Lya profile with a flux greater than an up- 
per limit value at that wavelength. Given the various 
combinations of Ig, FWHM, v outl N out , Th 21 and Nh 2 , 
a variety of Lya profiles adequately fit the Zj nc values. 
This is demonstrated in Figure [T] where several combi- 
nations of I Lya, I a b, and I inc are shown. We quantify 
the overall distribution of Ii, va profiles for each target 
by integrating each lL ya profile from 1210 A to 1222 A, 
resulting in fluxes Fz, ya (in units of ergs cm~ 2 s _1 ) which 
form an approximately Gaussian distribution. 

The average F^ya (< F^ya >) of each target is listed 
in Table [I] along with the la width of the distribution. 
< FLya > represents the total Lya flux at the star, 
and is generally constrained to within 10 - 20 %. We 
also include the average emission FWHM and H I ab- 
sorber properties (N out and v out ). Integrating I a b, the 
H2-incident Lya flux, yields < F a i, >. Assuming a ther- 
mally populated ground state distribution, typical Tjj 2 



values are -2500 ± 1000 K and log(A ^ (cm~ 2 )) -19 



± 1, similar to those found in TW Hya (Herczeg et al 
2004 1 . Strong UV and X-ray irradiation from the central 



star would act to preferentially populate excited rovibra- 
tion al levels in exc e ss of the thermal distribution (see 
g., Nomura et al. (2007[)). In our simple model, non- 



thermal excitation would be approximated by higher ro 
tational temperatures, and this increase in rovibrational 
population will be offset by lower total H2 column densi- 
ties. The varying effects of non-thermal excitation most 
likely contribute to the spread of the observed Nh 2 and 
Th 2 values. Lower Njj 2 can require increased Lya flux, 
although we expect that this increase will be of the same 
order of the uncertainty on the derived Lya emission flux 
(~ 10 - 20 %). Untangling the effects of non-thermal 
populations requires more complex modeling, including 
the effects of H2 formation on dust grains, and disk mod- 
eling efforts that predict the UV fluorescent spectrum of 
H 2 (and CO) would be very valuable. Regardless, the 
overall model Lya properties follow approximately gaus- 
sian distributions. Figure [2] compares < FL ya > and the 
observed Lya fluxes (i^o&s, summed over the entire pro- 
file) with the total H2 fluxes. F b s has a Spearman p 
rank correlation coefficient of 0.57 while < Fl vci >, by 
design, has a p of 0.96. 

Ideally we would compare our model Lya profiles with 
the observed Lya profiles. However, this involves ad- 
ditional parameters and a more complex analysis than 
presented here. Primarily, the radiative transfer of Lya 
photons between the source, the H2, and Earth must ac- 
count for the geometric filling fraction 77. This fraction 
represents the solid angle of H 2 illuminated by the Lya 
emission, such that r\=\ implies the H2 completely sur- 
rounds the star. An rj other than 1 would enter Equation 
2 by multiplying Ji„ c (A m ) by -q. Herczeg et al. (|2004| 
found 77=0.25 for T W Hya, imply i ng ga s in a dis^ 



k sur- 



face layer. However, Wood et al. (20021 measured r\ > 
2 in Mira B. A filling factor greater than 1 would possi- 
bly correspond to preferential scattering of H 2 photons 
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Figure 1. Observed and model spectra for V4046 Sgr (top row), RU Lupi (middle row), and DE Tau (bottom row). H2 emission lines 
are present throughout the FUV (left column), regardless of whether or not the photo-exciting Lya is observed (middle). In RU Lupi, a 
circumstellar wind extincts most of the short wavelength side of the Lya profile, while in DE Tau Lya is completely absorbed. We use 
the photo-excited H2 lines to reconstruct model Lya profiles (right). Each pair of intrinsic (I Lya, red) and outflow-absorbed (I a b, green) 
profiles are a best fit to a different set of [Nh 2 , Tf/ 2 ]-based incident fluxes (line, blue asterisks). 



into, or more likely Lya photons out of, our line of sight 



(Wood et al. 2002). Future analysis will incorporate r\ 
into a more detailed radiative transfer model of our sam- 
ple targets. 

3.3. Total Lya Flux in CTTSs 

We determine the fraction of the total FUV output 
from CTTSs that is in the Lya emission line, adopting a 
definition of the Lya fraction, fLya, given by 



fLya — 



F, 



Lya 



(F 



FUV 



''Lya) 



(5) 



where Fi, ya is our average model Lya flux from the star 
(in units of ergs cm~ 2 s _1 ; S3.2) and F FUV is the total 
unreddened FUV flux, measured over the 1160 - 1695 A 
bandpass (in units of ergs cm~ 2 s _1 ). The wings of the 
stellar Lya can extend to > ± 1000 km s _1 from line 
center in some targets (see e.g, Figure 2 (top) from Yang 
et al. 2011), therefore we excluded the 1207 - 1225 A 
region from the computation of Fpuy We also excluded 
the 1300 - 1310 A region to remove contamination by 
the geocoronal O I multiplet. A 5% error is assumed on 



Ffuv corresponding to the uncertainty in flux calibra- 
tion. The restriction of the total FUV band to 1160 - 
1695 A is mainly driven by the the lack of supporting 
FUSE observations at short wavelengths (912 - 1160 A) 
for all targets and the end of the STIS E140M bandpass 
at 1700 A. 

Molecular line e mission from H2 (Herczeg et al . 2006; 



this w ork) and CO (France et al. 2011a Schindhelm et al. 
2012 1 contributes to £'fuv- The FUV molecular emis- 



sion is thought to ori ginate from an inner warm disk 
surface (a < 3 AU) ([H erczeg ct al.||2004| |France et al 



2012[ ) or outflow ( |Saucedo et al |2003[| Walter et al.| 2003 
The total FUV flux also includes a molecular continuum 
whose excitation mech anism and spatial d i stribution are 
not we ll constrained (Bergin et al.||2004| |France et al. 



2011b I. < F^ya > shows a much better correlation with 
Fpuy than F b s (Figure [2]) . This is consistent with the 
idea that both the Lya ana FUV continuum emission are 
produced by related processes. 
We also calculate an observed Ly effraction (fobs) us- 

6s 



ing Fobs instead of Fl vci in Equation [5| We compare f a 



and fLya with the total H2 flux inTlgure [3] For our 
general assumption of r)=l, we find an average < fLya > 
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Figure 2. Integrated model (squares) and observed (triangles) Lya flux vs. total H2 flux (left) and total FUV flux (right) for each target. 
The total H2 flux is calculated by summing all the progression fluxes measured in our analysis. Fpjjy is the total unreddened FUV flux 
measured over the I f60 - 1695 A bandpass. 



of 81±9%, compared with for < f obs >. For rj = 

0.1, 0.25, 0.5, and 2.0, < f Lya > is 97%, 94%, 90%, and 
69%, respectively. A slight trend for decreasing < fhya > 
with increasing Fpuy is apparent, although this may be 
inappropriately weighted by BP Tau and TW Hya. Re- 
gardless, our calculated < fL ya > values demonstrate the 
validity of previous assertions that Lya emissio n domi 



nates the FU V radiation field from CTTSs (e.g. Bergin 
eFaL](p003|). 



4. DISCUSSION AND CONCLUSIONS 

We have demonstrated that Lya dominates the FUV 
emission of CTTSs. Lya- fluoresced H2 emission lines 
appear throughout the spectra of each CTTS in our sur- 
vey. Circumstellar and interstellar H I partially attenu- 
ates the line centers of the Lya profiles in most targets 
and completely absorbs them in several targets. The 
observed Lya fluxes do not correlate with either the to- 
tal H2 emission or the summed FUV continuum flux. 
The strongest fluorescent H2 progressions are used to re- 
construct the Lya profile incident on the molecular disk, 
yielding F^ ya values accurate to within ~ 20% in most 
targets. This uncertainty is caused by the distribution of 
possible H I and H2 properties. Our model assumptions 
(such as a single Gaussian emission component or sim- 
ple outflow absorber) may hide larger systematic errors, 
however future work will study these effects in a more de- 
tailed model. We find that the intrinsic Lya comprises 
81±9% of the total FUV emission from CTTSs, com- 
pared with a fraction of only 15^ 15 % for the observed 
Lya profiles. This demonstrates the need for Lya recon- 
struction to achieve accurate disk models of CTTSs. 



It is clear from our results that the detection of the 
strong Lya line in TW Hya was not limited to a single 
object. Rather, strong Lya emission dominates the FUV 
flux from all accreting (classical) T Tauri stars. Our 
measurements of the relative Lya/FUV continuum flux 
only compare Lya to the FUV flux from 1160 - 1695 A. 
Thus if there is significant UV flux shortward of 1160 A 
then the strength of Lya relative to the FUV radiation 
could go down. However, in TW Hya the flux below 
this limit is only 5% of the FUV flux below 1700 A. 
Similarly, France et al. (2011b) have shown that the FUV 
continuum decreases to shorter wavelengths across the 
FUV bandpass. In addition, the absorption properties 
of grains strongly limit the propagation of UV photons 
near the Lyman limit. 

The derived Lya fractions confirm the dominance of 
Lya in the FUV spectrum of the accreting young stars 
with disks. This is important because Lya photons from 
the star will see the atomic hydrogen layer, with shallow 
angle of incidence, above the molecular surface. Isotropic 
scattering will lead a significant fraction of the Lya flux 
to rain down on the disk with greater penetrating power 
than typical UV continuum photons ( Bethell & Bergin 



2011 1. In addition, the Lya emission reprocessed by rl 2 
(which is scattered throughout the FUV spectrum) wil l 



also emit directly on the disk surface ( France et al.|2012 ) 
These two effects increase the penetration ot U V pho 
tons beyond the simple case where one assumes the UV 
photons observe the disk surface with shallow angle of 
incidence and with the propagation solely influenced by 
grains. In general this should lead to greater heating and 
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Total H 2 flux (ergs s 



Figure 3. Observed (top) and model (bottom) Lycr fraction (as- 
suming r)=l) vs. total H2 flux. The diamond in each plot desig- 
nates the average fL ya value for the observed and model profiles. 
The significant increase in Lya fraction from the observed to the 
model profiles demonstrates the importance of Lyo reconstruction. 



additional chemical effects deeper in the disk system. 
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